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Abstract 

ff^ We review current progress in the functional renormalization group treatment of disordered systems. After an el- 

O ementary introduction into the phenomenology, we show why in the context of disordered systems a functional 

■ renormalization group treatment is necessary, contrary to pure systems, where renormalization of a single coupling 

^ constant is sufficient. This leads to a disorder distribution, which after a finite renomalization becomes non-analytic, 

(]J thus overcoming the predictions of the seemingly exact dimensional reduction. We discuss, how a renormalizable 

field theory can be constructed, even beyond 1-loop order. We then discuss an elastic manifold imbedded in N di- 

mensions, and give the exact solution for ^ cxd. This is compared to predictions of the Gaussian replica variational 

^ ansatz, using replica symmetry breaking. We finally discuss depinning, both isotropic and anisotropic, and the scaling 

^ function for the width distribution of an interface. 
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1 Introduction 



Statistical mechanics is by now a rather mature branch of physics. For pure systems Uke a ferromagnet, it allows 
to calculate so precise details as the behavior of the specific heat on approaching the Curie-point. We know that it 
diverges as a function of the distance in temperature to the Curie-temperature, we know that this divergence has the 
form of a power-law, we can calculate the exponent, and we can do this with at least 3 digits of accuracy. Best of all, 
these findings are in excellent agreement with the most precise experiments. This is a true success story of statistical 
mechanics. On the other hand, in nature no system is really pure, i.e. without at least some disorder ("dirt"). As 
experiments (and theory) seem to suggest, a little bit of disorder does not change the behavior much. Otherwise 
experiments on the specific heat of Helium would not so extraordinarily well confirm theoretical predictions. But 
what happens for strong disorder? By this I mean that disorder completely dominates over entropy. Then already the 
question: "What is the ground-state?" is no longer simple. This goes hand in hand with the appearance of so-called 
metastable states. States, which in energy are very close to the ground-state, but which in configuration-space may 
be far apart. Any relaxational dynamics will take an enormous time to find the correct ground-state, and may fail 
altogether, as can be seen in computer-simulations as well as in experiments. This means that our way of thinking, 
taught in the treatment of pure systems, has to be adapted to account for disorder. We will see that in contrast to 
pure systems, whose universal large-scale properties can be described by very few parameters, disordered systems 
demand the knowledge of the whole disorder-distribution function (in contrast to its first few moments). We show 
how universality nevertheless emerges. 

Experimental realizations of strongly disordered systems are glasses, or more specifically spin-glasses, vortex- 
glasses, electron-glasses and structural glasses (not treated here). Furthermore random-field magnets, and last not 
least elastic systems in disorder. 

What is our current understanding of disordered systems? It is here that the success story of statistical mechanics, 
with which I started, comes to an end: Despite 30 years of research, we do not know much: There are a few exact 
solutions, there are phenomenological methods (like the droplet-model), and there is the mean-field approximation, 
involving a method called replica-symmetry breaking (RSB). This method is correct for infinitely connected systems, 
e.g. the SK-model (Sherrington Kirkpatrick model), or for systems with infinitely many components. However it is 
unclear, to which extend it applies to real physical systems, in which each degree of freedom is only coupled to a 
finite number of other degrees of freedom. 

Another interesting system are elastic manifolds in a random media, which has the advantage of being approach- 
able by other (analytic) methods, while still retaining all the rich physics of strongly disordered systems. Here, I 
review recent advances obtained in collaboration with Pierre Le Doussal [1-12]. This review is an extended version 
of [13]. 

2 Physical realizations, model and observables 

Before developing the theory to treat elastic systems in a disordered environment, let us give some physical realiza- 
tions. The simplest one is an Ising magnet. Imposing boundary conditions with all spins up at the upper and all spins 
down at the lower boundary (see figure 1), at low temperatures, a domain wall separates a region with spin up from 
a region with spin down. In a pure system at temperature T = 0, this domain wall is completely flat. Disorder can 
deform the domain wall, making it eventually rough again. Two types of disorder are common: random bond (which 
on a course-grained level also represents missing spins) and random field (coupling of the spins to an external random 
magnetic field). Figure 1 shows, how the domain wall is described by a displacement field u{x). Another example 
is the contact line of water (or liquid Helium), wetting a rough substrate, see figure 2. (The elasticity is long range). 
A realization with a 2-parameter displacement field u{x) is the deformation of a vortex lattice: the position of each 
vortex is deformed from x to x + u{x). A 3-dimensional example are charge density waves. 
All these models have in common, that they are described by a displacement field 

xeW^ — ^ u{x) e . (2.1) 
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Figure 1: An Ising magnet at low temperatures forms a domain wall described by a function u{x) (right). An experiment on a thin Cobalt film 
(left) [14]; with kind permission of the authors. 



For simplicity, we set = 1 in the following. After some initial coarse-graining, the energy ?Y — + Wdo 
consists out of two parts: the elastic energy 



and the disorder 



HdoM = / d'^xVix,uix)) . 



(2.2) 



(2.3) 



In order to proceed, we need to specify the correlations of disorder. Suppose that fluctuations u in the transversal 
direction scale as 

{u{x) - u{y)f \x - y\^'^ (2.4) 

with a roughness-exponent C < 1- Starting from a disorder correlator 



V{u, x)V{u', x') = f{x - x')R{u - u') 



(2.5) 




Figure 2: A contact line for the wetting of a disordered substrate by Glycerine [15]. Experimental setup (left). The disorder consists of 
randomly deposited islands of Chromium, appearing as bright spots (top right). Temporal evolution of the retreating contact-line (bottom 
right). Note the different scales parallel and perpendicular to the contact-line. Pictures courtesy of S. Moulinet, with kind permission. 
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Figure 3: Cartoon of an elastic lattice (e.g. vortex lattice) deformed by disorder. This is described by a vector u{x). 

and performing one step in the RG-procedure, one has to rescale more in the x-direction than in the u-direction. This 
will eventually reduce f{x — x') to a ^-distribution, whereas the structure of R{u — u') remains visible. We therefore 
choose as our starting-model 

V{u,x)V{u',x') := 6'^{x - x')R{u - u') . (2.6) 

There are a couple of useful observables. We already mentioned the roughness-exponent C,. The second is the 
renormalized (effective) disorder. It will turn out that we actually have to keep the whole disorder distribution function 
R{u), in contrast to keeping a few moments. Other observables are higher correlation functions or the free energy. 

3 Treatment of disorder 

Having defined our model, we can now turn to the treatment of disorder. The problem is to average not the partition- 
function, but the free energy over disorder: = In Z. This can be achieved by the beautiful replica-trick. The idea 
is to write 

InZ = lim - fe"'""^ - O = Hm - (Z" - 1) (3.1) 

n^O n \ / n^O n 

and to interpret as the partition-function of an n times replicated system. Averaging e~^'»=i^° over disorder 
then leads to the replica-Hamiltonian 

^N = ^X; / d"xi(Vn,(x))2-^ ^ I d''xRMx)-u,{x)). (3.2) 

a=l a,b=l 

Let us stress that one could equivalently pursue a dynamic or a supersymmetric formulation. We therefore should 
not, and in fact do not encounter, problems associated with the use of the replica-trick. 

4 Dimensional reduction 

There is a beautiful and rather mind-boggling theorem relating disordered systems to pure systems (i.e. without 
disorder), which applies to a large class of systems, e.g. random field systems and elastic manifolds in disorder. It is 
called dimensional reduction and reads as follows [16]: 

Theorem: A d-dimensional disordered system at zero temperature is equivalent to all orders in perturbation theory 
to a pure system in d — 2 dimensions at finite temperature. Moreover the temperature is (up to a constant) nothing 
but the width of the disorder distribution. A simple example is the 3-dimensional random-field Ising model at zero 
temperature; according to the theorem it should be equivalent to the pure 1 -dimensional Ising-model at finite temper- 
ature. But it has been shown rigorously, that the former has an ordered phase, whereas we have all solved the latter 
and we know that there is no such phase at finite temperature. So what went wrong? Let me stress that there are no 
missing diagrams or any such thing, but that the problem is more fundamental: As we will see later, the proof makes 
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Figure 4: The disorder vertex R{ua{x) — Ub{x)) and the correlation-function C{x — y), with Fourier-transform C{k) = which is 
diagonal in replica-space (left). Contracting two disorder-vertices with two correlation-functions leads to the two 1-loop contributions SR to 
the disorder-correlator R (right). The integral C{x — yf = where L is some IR-cutoff. 

assumptions, which are not satisfied. Nevertheless, the above theorem remains important since it has a devastating 
consequence for all perturbative calculations in the disorder: However clever a procedure we invent, as long as we do 
a perturbative expansion, expanding the disorder in its moments, all our efforts are futile: dimensional reduction tells 
us that we get a trivial and unphysical result. Before we try to understand why this is so and how to overcome it, let 
me give one more example. Dimensional reduction allows to calculate the roughness-exponent ( defined in equation 
(2.4). We know (this can be inferred from power-counting) that the width n of a d-dimensional manifold at finite 
temperature in the absence of disorder scales as u ~ x^'^~^^^'^. Making the dimensional shift implied by dimensional 
reduction leads to 

{u{x) - u{0)f ~ x^-*^ = i.e. C = ^^- (4.1) 

5 The Larkin-length 

To understand the failure of dimensional reduction, let us turn to an interesting argument given by Larkin [17]. He 
considers a piece of an elastic manifold of size L. If the disorder has correlation length r, and characteristic potential 
energy /, this piece will typically see a potential energy of strength 

d 

^Do=/(^y . (5.1) 

On the other hand, there is an elastic energy, which scales like 



Eei = cL'^-2. (5.2) 



These energies are balanced at the Larkin-length L = Lc with 



2 \ 4-d 



i^c = ( ) . (5.3) 



More important than this value is the observation that in all physically interesting dimensions d < 4, and at scales 
L > Lc, the membrane is pinned by disorder; whereas on small scales elastic energy dominates. Since the disorder 
has a lot of minima which are far apart in configurational space but close in energy (metastabiUty), the manifold can 
be in either of these minimas, and the ground-state is no longer unique. However exactly this is assumed in e.g. the 
proof of dimensional reduction; as is most easily seen in its supersymmettic formulation [18]. 

6 The functional renormalization group (FRG) 

Let us now discuss a way out of the dilemma: Larkin's argument suggests that 4 is the upper critical dimension. So 
we would like to make an e = 4 — d expansion. On the other hand, dimensional reduction tells us that the roughness is 
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( = (see (4.1)). Even though this is systematically wrong below four dimensions, it tells us correctly that at the 
critical dimension d = A, where disorder is marginally relevant, the field u is drmensionless. This means that having 
identified any relevant or marginal perturbation (as the disorder), we find immediately another such perturbation by 
adding more powers of the field. We can thus not restrict ourselves to keeping solely the first moments of the disorder, 
but have to keep the whole disorder-distribution function R(u). Thus we need afunctional renormalization group 
treatment (FRG). Functional renormalization is an old idea going back to the seventies, and can e.g. be found in [19]. 
For disordered systems, it was first proposed in 1986 by D. Fisher [20]. Performing an infinitesimal renormalization, 
i.e. integrating over a momentum shell a la Wilson, leads to the flow dtR{u), with (e = 4 — d) 

diR{u) = (e - 4C) R{u) + CuR'{u) + ^R"{uf - R"{u)R"{0) . (6.1) 

The first two terms come from the rescaUng of R and u respectively. The last two terms are the result of the 1-loop 
calculations, which are sketched in figure 4. 

More important than the form of this equation is it actual solution, sketched in figure 5. After some finite renor- 
malization, the second derivative of the disorder R"{u) acquires a cusp at n = 0; the length at which this happens is 
the Larkin-length. How does this overcome dimensional reduction? To understand this, it is interesting to study the 
flow of the second and forth moment. Taking derivatives of (6.1) w.r.t. u and setting u to 0, we obtain 

5^i?"(0) = (e - 20 i2"(0) + i2'"(0)2 (e-2C)i?"(0) (6.2) 
diR!"'{0) = eR""{0) + 3R""{0f + 4R"'{0)R""'{0) eR""{0) + 3R""{0f . (6.3) 

Since R{u) is an even function, R"'{0) and R""'{0) are 0, which we have already indicated in Eqs. (6.2) and (6.3) . 
The above equations for R"{0) and R""{0) are in fact closed. Equation (6.2) tells us that the flow of R"{0) is trivial 
and that ( = e/2 = This is exactly the result predicted by dimensional reduction. The appearance of the cusp 
can be inferred from equation (6.3). Its solution is 

Thus after a finite renormalization R""{0) becomes infinite: The cusp appears. By analyzing the solution of the 
flow-equation (6.1), one also finds that beyond the Larkin-length R"{0) is no longer given by (6.2) with i?'"(0)^ = 0. 
The correct interpretation of (6.2), which remains vaUd after the cusp-formation, is (for details see below) 

deR"{0) = (e - 20 R"{0) + R"'{0+f . (6.5) 

Renormalization of the whole function thus overcomes dimensional reduction. The appearance of the cusp also 
explains why dimensional reduction breaks down. The simplest way to see this is by redoing the proof for elastic 
manifolds in disorder, which in the absence of disorder is a simple Gaussian theory. Terms contributing to the 2-point 
function involve R"{0), TR""{0) and higher derivatives of R{u) atu = 0, which all come with higher powers of T. 
To obtain the limit of T ^ 0, one sets T = 0, and only R"{0) remains. This is the dimensional reduction result. 
However we just saw that R""(0) becomes infinite. Thus R""{0)T may also contribute, and the proof fails. 
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Figure 6; The toy model (7.1) before (left) and after (right) the Larkin-scale. 



7 Why is a cusp necessary? 

The appearance of a cusp might suggest that this approach is fatally ill. Let me present a simple argument due to 
Leon Balents [21], why a cusp is a physical necessity and not an artifact. To this aim, consider a toy model with only 
one Fourier-mode u = Ug 

n[u] = ^q''u^ + V'eV{u). (7.1) 

Since equation (6.1) has a fixed point of order R{u) ~ e for all e > 0, V{u) scales like ^/e for e small and we have 
made this dependence explicit in (7.1) by using V{u) = ^/eV{u). The only further input comes from the physics: 
For L < Lc, i.e. before we reach the Larkin length, there is only one minimum, as depicted in figure 6. On the other 
hand, for L > Lc, there are several minima. Thus there is at least one point for which 

^ n[u] =q^ + ^feViu) < . (7.2) 

In the limit of e — > 0, this is possible if and only if -R""{0), which a priori should be finite for e — > 0, becomes 
infinite: 

-^R""{0) = V"{u)V"{u')l^^, = oo . (7.3) 
This argument shows that a cusp is indeed a physical necessity. 



8 Beyond 1 loop? 

Functional renormalization has successfully been applied to a bunch of problems at 1-loop order. From a field theory, 
we however demand more. Namely that it 

• allows for systematic corrections beyond 1-loop order 

• be renormalizable 

• and thus allows to make universal predictions. 

However, this has been a puzzle since 1986, and it has even been suggested that the theory is not renormalizable due 

3 

to the appearance of terms of order £2 [22]. Why is the next order so complicated? The reason is that it involves 
terms proportional to R"'{0). A look at figure 3 explains the puzzle. Shall we use the symmetry of R{u) to conclude 
that R"'{0) is 0? Or shall we take the left-hand or right-hand derivatives, related by 

i?'"(0+) := lim R"'{u) = - lim R"'{u) =: -R"'{0-). (8.1) 

In the following, I will present our solution of this puzzle, obtained at 2-loop order in [1] and then at large in [2]. 
The latter approach allows for another independent control-parameter, and sheds further light on the cusp-formation. 
We refer the reader to these works for further references. 
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Figure 7: Results for C, in the random bond case. 

9 Results at 2-loop order 

For the flow-equation at 2-loop order, we find [1] 

diR{u) = (e - 40 R{u) + QuR'{u) + ]^R"{uf - R"{u)R"{0) 
+i _ i?"(o)) R'"(uf - ^R"'{0+fR"iu) . 



(9.1) 



The first line is the result at 1-loop order, already given in (6.1). The second line is new. The most interesting 
term is the last one, which involves i?'"(0"'")^ and which we therefore call anomalous. The hard task is to fix the 
prefactor (— ^). We have up to now invented six algorithms to do so; one leads to inconsistencies and shall not be 
reported here. The other five algorithms are consistent with each other: The sloop-algorithm, recursive construction, 
reparametrization invariance, renormalizability, and potentiality [1,4]. For lack of space, we restrain our discussion 
to the last two ones. At 2-loop order the following diagram appears 




R"{0)) R!"{uf 



]^R:'{u)R"'{Q+f 



(9.2) 



leading to the anomalous term. The integral (not written here) contains a subdivergence, which is indicated by the 
box. Renormalizability demands that its leading divergence (which is of order 1/e^) be canceled by a 1-loop counter- 
term. The latter is unique thus fixing the prefactor of the anomalous term. (The idea is to take the 1-loop correction 
5R in figure 2 and replace one of the R" in it by 5R" itself, which the reader can check to leading to the terms given 
in (9.2) plus terms which only involve even derivatives.) 

Another very physical demand is that the problem remain potential, i.e. that forces still derive from a potential. 
The force-force correlation function being —R"{u), this means that the flow of ii'(O) has to be strictly 0. (The 
simplest way to see this is to study a periodic potential.) From (9) one can check that this does not remain true if one 
changes the prefactor of the last term in (9); thus fixing it. 

Let us give some results for cases of physical interest. First of all, in the case of a periodic potential, which is 
relevant for charge-density waves, the fixed-point function can be calculated analytically as (we choose period 1, the 
following is for u G [0, 1]) 



R*{u) 



e 

72 



uHi 



const. 



(9.3) 



This leads to a universal amplitude. In the case of random field disorder (short-ranged force-force correlation func- 
tion) C = §■ For random-bond disorder (short-ranged potential-potential correlation function) we have to solve (9.1) 
numerically, with the result C, = 0.20829804e + 0.006858e^. This compares well with numerical simulations, see 
figure 7. 



10 Large 

In the last section, we have discussed renormalization in a loop expansion, i.e. expansion in e. In order to indepen- 
dently check consistency it is good to have a non-perturbative approach. This is achieved by the large- limit, which 
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can be solved analytically and to which we turn now. We start from 



1 " /• 

a,fe=l ^ 

where in contrast to (3.2), we use an A'^-component field u. For = 1, we identify B{u^) = R{u). We also have 
added a mass m to regularize the theory in the infra-red and a source j to calculate the effective action T{u) via a 
Legendre transform. For large AT the saddle point equation reads [2] 

B'{ul,) = B' (ul^ + 2Th + 4/2[5'(u2,) - ^'(0)]) (10.2) 

This equation gives the derivative of the effective (renormalized) disorder S as a function of the (constant) background 
field -u^^ = {ua — Ub)'^ in terms of: the derivative of the microscopic (bare) disorder B, the temperature T and the 

integrals /„ := (p^^. 

The saddle-point equation can again be turned into a closed functional renormahzation group equation for B by 
taking the derivative w.r.t. m: 

mr) 1 fTR'(T\ 

diB{x) ^ -Q^B{x) = (e - 4C)Bix) + 2^xB'{x) + -B'{xf - B'{x)B'{0) + ^^^^ (10.3) 

This is a complicated nonlinear partial differential equation. It is therefore surprising, that one can find an analytic 
solution. (The trick is to write down the flow-equation for the inverse function of B\x), which is linear.) Let us only 
give the results of this analytic solution: First of all, for long-range correlated disorder of the form B'{x) ~ x'"^ , the 
exponent C can be calculated analytically as C = 2(1+^) • It agrees with the replica-treatment in [24] and the 1-loop 
treatment in [22]. For short-range correlated disorder, C = 0. Second, it demonstrates that before the Larkin-length, 
B{x) is analytic and thus dimensional reduction holds. Beyond the Larkin length, B"{0) = 00, a cusp appears and 
dimensional reduction is incorrect. This shows again that the cusp is not an artifact of the perturbative expansion, but 
an important property even of the exact solution of the problem (here in the Umit of large N). 



1 1 Relation to Replica Symmetry Breaking (RSB) 

There is another treatment of the limit of large given by Mezard and Parisi [24]. They start from (10.1) but without 
a source-term j. In the Umit of large N, a Gaussian variational ansatz of the form 

a=l •'^ 0,6=1 

becomes exact. The art is to make an appropriate ansatz for a ah- The simplest possibility, aab = f for all a 7^ 6 
reproduces the dimensional reduction result, which breaks down at the Larkin length. Beyond that scale, a replica 
symmetry broken (RSB) ansatz for Gah is suggestive. To this aim, one can break Gah into four blocks of equal size, 
choose one (variationally optimized) value for the both outer diagonal blocks, and then iterate the procedure on the 
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diagonal blocks, resulting in 



\ 



<7ab = 



\ 




I 



(11.2) 



One finds that the more often one iterates, the better the results. In fact, one has to repeat this procedure infinite 
many times. This seems like a hopeless endeavor, but Parisi has shown that the infinitely often replica symmetry 
broken matrix can be parameterized by a function [o"](z) with z G [0, 1]. In the SK-model, z has the interpretation of 
an overlap between replicas. While there is no such simple interpretation for the model (11.1), we retain that z = 
describes distant states, whereas z = 1 describes nearby states. The solution of the large- saddle-point equations 
leads to the curve depicted in figure 6. Knowing it, the 2-point function is given by 



The important question is: What is the relation between the two approaches, which both pretend to calculate the same 
2-point function? Comparing the analytical solutions, we find that the 2-point function given by FRG is the same as 
that of RSB, if in the latter expression we only take into account the contribution from the most distant states, i.e. those 
for z between and z^ (see figure 8). To understand why this is so, we have to remember that the two calculations 
were done under quite different assumptions: In contrast to the RSB-calculation, the FRG-approach calculated the 
partition function in presence of an external field j, which was then used to give via a Legendre transformation the 
effective action. Even if the field j is finally turned to 0, the system might remember its preparation, as is the case 
for a magnet: Preparing the system in presence of a magnetic field will result in a magnetization which aligns with 
this field. The magnetization will remain, even if finally the field is turned off. The same phenomena happens here: 
By explicitly breaking the replica-symmetry through an applied field, all replicas will settle in distant states, and the 
close states from the Parisi-function [a] (z) + (which describes spontaneous RSB) will not contribute. However, 
we found that the full RSB -result can be reconstructed by remarking that the part of the curve between Zm and Zc is 
independent of the infrared cutoff m, and then integrating over m [2] (rric is the mass corresponding to Zc): 




(11.3) 




(11.4) 



[a]{z) + m 



from UV-cutoff 



IR-cutoff 





2 




m 







z 



m 



1 



Figure 8: The function [a] (u) + as given in [24]. 
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We also note that a similar effective action has been proposed in [25]. While it agrees qualitatively, it does not 
reproduce the correct FRG 2-point function, as it should. 



12 Finite 



Up to now, we have studied the functional RG in two cases: For one component = 1 and in the limit of a large 
number of components, ^ oo. The general case of finite is more difficult to handle, since derivatives of the 
renormalized disorder now depend on the direction, in which this derivative are taken. Define amplitude u := \u\ 
and direction u := u/\u\ of the field. Then deriving the latter variable leads to terms proportional to 1/u, which are 
diverging in the limit of u ^ 0. This poses additional problems in the calculation, and it is a priori not clear that the 
theory at A^ / 1 exists, supposed this is the case for A^ = 1. At 1-loop order everything is well-defined [22]. We 
have found a consistent RG-equation at 2-loop order [8] : 



The first line is the 1-loop equation, given in [22]. The second and third line represent the 2-loop equation, with the 
new anomalous terms proportional to i?'"(0+)^ (third line). 

The fixed point equation (12.1) can be integrated numerically, order by order in e. The result, specialized to 
directed polymers, i.e. e = 3 is plotted on figure 9. We see that the 2-loop corrections are rather big at large A^, so 
some doubt on the applicability of the latter down to e = 3 is advised. However both 1- and 2-loop results reproduce 
well the two known points on the curve: C = 2/3forA^=l and C = for A^ = oo. The latter result has been 
given in section 10. Via the equivalence [26] of the directed polymer problem in A^ dimensions treated here and the 
KPZ-equation of non-linear surface growth in A^ dimensions, which relate the roughness exponent C of the directed 
polymer to the dynamic exponent zrpz in the KPZ-equation via ( = j^^, we know that ({N = 1) = 2/3. 

The line ( = 1/2 given on figure 9 plays a special role: In the presence of thermal fluctuations, we expect the 
roughness-exponent of the directed polymer to be bounded hy ( > 1/2. In the KPZ-equation, this corresponds to 
a dynamic exponent zrpz = 2, which via the exact scaling relation zrpz + Ckpz = 2 is an upper bound in the 
strong-coupling phase. The above data thus strongly suggest that there exists an upper critical dimension in the KPZ- 
problem, with due ~ 2.4. Even though the latter value might be an underestimation, it is hard to imagine what can 
go wrong qualitatively with this scenario. The strongest objections will probably arize from numerical simulations, 





2 N-1 {R'{u)-uR"{u)f {2R'{u)+u{R"iu)-3R"{0))) 
N -lR'{ry 



(12.1) 
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Figure 9: Results for the roughness <^ at 1- and 2-loop order, as a function of the number of components A'^. 
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such as [27]. However the latter use a discrete RSOS model, and the exponents are measured for interfaces, which in 
large dimensions have the thickness of the discretization size, suggestions that the data are far from the asymptotic 
regime. We thus strongly incourage better numerical simulations on a continuous model, in order to settle this issue. 

13 Depinning transition 

Another important class of phenomena for elastic manifolds in disorder is the so-called "depinning transition": Ap- 
plying a constant force to the elastic manifold, e.g. a constant magnetic field to the ferromagnet mentioned in the 
introduction, the latter will only move, if a certain critical threshhold force /c is surpassed, see figure 10. (This is 
fortunate, since otherwise the magnetic domain walls in the hard-disc drive onto which this article is stored would 
move with the effect of deleting all information, depriving the reader from his reading.) At / = fc, the so-called 
depinning transition, the manifold has a distinctly different roughness exponent C (see Eq. (2.4)) from the equilibrium 
(/ = 0). For / > fc, the manifold moves, and close to the transition, new observables and corresponding exponents 
appear: 

• the dynamic exponent z relating correlation functions in spatial and temporal direction 

• a correlation length set by the distance to fc 

^-\f-fcr 

• furthermore, the new exponents are not all independent, but satisfy the following exponent relations [28] 

P = u{z-C) ^=2^ ^^^-^^ 

The equation describing the movement of the interface is 

dtu{x,t) = (V2 + m'^)u{x,t) + F{u{x,t), x) , F{u,x) = —duV{u,x) (13.2) 

This model has been treated at 1-loop order by Natterman et al. [28] and by Narayan and Fisher [29]. The 1-loop 
flow-equations aie identical to those of the statics. This is surprising, since physically, the phenomena at equilibrium 
and at depinning are quite different. There is even the claim by [29], that the roughness exponent in the random 
field universality class is exactly C = e/3, as it is in the equilibrium random field class. After a long debate among 
numerical physicists, the issue is today resolved: The roughness is significantly larger, and reads e.g. for the driven 
polymer = 1.25, instead of = 1 as predicted in [29]. Clearly, a 2-loop analysis [3] is necessary, to resolve these 
issues. Such a treatment starts from the dynamic action 




(13.3) 



fc f 

Figure 10: Velocity of a pinned interface as a function of the applied force. Zero force: equilibrium, f = fc'- depinning. 
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where the "response field" u{x,t) enforces the equation of motion (13.2) and 

F{u, x)F{u',x') = A{u - u')S'^{x - x') = -R"{u - u')5'^{x - x') (13.4) 

is the force-force correlator, leading to the second term in (13.3). As in the statics, one encounters terms proportional 
to A'(0+) = — ii!'"(0+). Here the sign-problem can uniquely be resolved by supposing that the membrane only 
jumps ahead, 

t>t' u{x,t) > u{x,t') . (13.5) 

Practically this means that when evaluating diagrams containing A{u{x, t) — u{x, t')), one splits them into two pieces, 
one with t < t' and one with t > t'. Both pieces are well defined, even in the limit of t ^ t'. As the only tread-off of 
this method, diagrams can become complicated and difficult to evaluate; however they are always well-defined. 

Physically, this means that we approach the depinning transition from above. This is reflected in (13.5) by the 
fact that u{x, t) may remain constant; and indeed correlation-functions at the depinning transition are completely 
independent of time [3]. On the other hand a theory for the approach of the depinning transition from below (/ < /c) 
has been elusive so far. 

At the depinning transition, the 2-loop functional RG reads [1,3] 

deR{u) = {e - 4C)R{u) + CuR'{u) + ^R"{uf-R"{u)R"{0) 

+ i [R"{u) - R"{0)] R"'{uf + ^ R"'{id+fR"{u) (13.6) 

First of all, note that it is a priori not clear that the functional RG equation, which is a flow equation for A(u) = 
—R"{u) can be integrated to a functional RG-equation for R{u). We have chosen this representation here, in order 
to make the difference to the statics evident: The only change is in the last sign on the second line of (13.6). This has 
important consequences for the physics: First of all, the roughness exponent ( for the random-field universaUty class 
changes from ^ = | to 

( = l(l + 0.14331e-F...) (13.7) 

Second, the random-bond universality class is unstable and always renormalizes to the random-field universality class, 
as is physically expected: Since the membrane only jumps ahead, it always experiences a new disorder configuration, 
and there is no way to know of wether this disorder can be derived from a potential or not. This non-potentiality 
is most strikingly observed in the random periodic universaUty class, which is the relevannt one for charge density 
waves. The fixed point for a periodic disorder of period one reads (remember A(n) = —R"{u)) 

€ e2 /e £2 



^>) = 36^-108 - 6+9)^^^-^) ^''-'^ 





d 


e 




estimate 


simulation 




3 


0.33 


0.38 


0.38±0.02 


0.34±0.01 


c 


2 


0.67 


0.86 


0.82±0.1 


0.75±0.02 




1 


1.00 


1.43 


1.2±0.2 


1.25±0.01 




3 


0.89 


0.85 


0.84±0.01 


0.84±0.02 


(3 


2 


0.78 


0.62 


0.53±0.15 


0.64±0.02 




1 


0.67 


0.31 


0.2±0.2 


0.25 ...0.4 




3 


0.58 


0.61 


0.62±0.01 




V 


2 


0.67 


0.77 


0.85±0.1 


0.77±0.04 




1 


0.75 


0.98 


1.25±0.3 


1±0.05 



Figure 11: The critical exponents at the depinning transition, 
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Integrating over a period, we find (suppressing in F{u, x) the dependence on the coordinate x for simplicity of 
notation) 

duA*{u) = / duF{u)F{u') = -j^ ■ (13.9) 
Jo 10° 

In an equiUbrium situation, this correlator would vanish, since potentiahty requires dn F{u) = 0. Here, there are 
non-trivial contributions at 2-loop order (order e^), violating this condition and rendering the problem non-potential. 
This same mechanism is also responsible for the violation of the conjecture C = f , which could be proven on the 
assumption that the problem remains potential under renormalization. Let us stress that the breaking of potentiality 
under renormalization is a quite novel observation here. 

The other critical exponents mentioned above can also be calculated. The dynamical exponent z (for RF-disorder) 
reads [1,3] 

z = 2-^e-0.04321e2 + ... (13.10) 

All other exponents are related via the relation (13.1). That the method works well even quantitatively can be infered 
from figure 1 1 . 



14 Universal width distribution 



Exponents are not the only interesting quantities: In experiments and simulations, often whole distributions can be 
measured, as e.g. the universal width distribution of an interface that we have computed at depinning [12,11]. Be (u) 
the average position of an interface for a given disorder configuration, then the spatially averaged width 



w 



j-^ Iju{x) - {u)y 



(14.1) 



is a random variable, and we can try to calculate and measure its distribution P{w'^). The rescaled function ^{z), 
defined by 

P{w'^) = 1/^$ (w'^/llp) (14.2) 



will be universal, i.e. independent of microscopic details and the size of the system. 

Supposing all correlations to be Gaussian, ^{z) can be calculated analytically. It depends on two parameters, 
the roughness exponent ^ and the dimension d. Numerical simulations displayed on figure 12 show spectacular 



analytic result 
simulation, L = 256 




Figure 12: Scaling function ^(z) for the (1 + l)-diniensional harmonic model, compared to the Gaussian approximation for = 1.25. Data 
from [12]. 
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agreement between analytical and numerical results. As expected, the Gaussian approximation is not exact, but to 
see deviations in a simulation, about 10^ samples have to be used. Analytically, corrections can be calculated: They 
are of order R"'{0^)^ and small. Physically, the distribution is narrower than a Gaussian. 

15 Anisotropic depinning, directed percolation, branching and all that 

We have discussed in section 13 isotopic depinning, which as the name suggest is a situation, where the system is 
invariant under a tilt. This isotropy can be broken through an additional anharmonic elasticity 

£^ela.tic = [Vu{x)f + C4 [Vu{x)f , (15.1) 

leading to a drastically different universality class, the so-called anisotropic depinning universaUty class, as found 
recently in numerical simulations [30]. It has been observed in simulations [31,32], that the drift- velocity of an 
interface is increased, which can be intepreted as a tilt-dependent term, leading to the equation of motion in the form 

dtu{x, t) = V'^u{x, t)+X [Vu{x, t)f + F{x, u{x, t)) + f . (15.2) 

However it was for a long time unclear, how this new term (proportional to A), which usually is refered to as a KPZ- 
term, is generated, especially in the limit of vanishing drift-velocity. In [5] we have shown that this is possible in a 
non-analytic theory, due to the diagram given in figure 13. 

For anisotropic depinning, numerical simulations based on cellular automaton models which are believed to be 
in the same universaUty class [33,34], indicate a roughness exponent ( 0.63 in d = 1 and ^ ~ 0.48 in d = 2. On 
a phenomenological level it has been argued [33-35] that configurations at depinning can be mapped onto directed 
percolation in d = 1 dimensions, which yields indeed a roughness exponent Cdp = ^±/i^\\ = 0.630 ± 0.001, and 
it would be intriguing to understand this from a systematic field theory. 

This theory was developed in [5], and we review the main resulsts here. A strong simpUfication is obtained by 
going to the Cole-Hopf transformed fields 

:=e^"(^'*) ^ = M^i^ . (15.3) 

A 

The equation of motion becomes after multiplying with XZ{x, t) (dropping the term proportional to /) 

dtZ{x, t) = V^Z{x, t) + XF H^i^i^ j z{x, t) (15.4) 

and the dynamical action (after averaging over disorder) 

S= j Z{x,t){dt-V^)Z{x,t)-% [ Z{x, t)Z{x, t) A f ~ ^'^ \ Z{x, t')Z{x, t') (15.5) 

Jxt 2 J^ff, \ X J 

This leads to the FRG flow equation at 1-loop order 

deA{u) = (e - 2QA{u) + CuA'{u) - A"{u) {A{u) - A(0)) - A'{uf 

+2AA(u)A'(0+) + 2A2 {A{uf + A(n) A(0)) (15.6) 



Figure 13: The diagram generating the irreversible nonlinear KPZ term with one disorder vertex and one C4 vertex (the bars denote spatial 
derivatives). 
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Figure 14: Fixed point structure for A = 2, which is a typical value. The ratio c/fe is not renormalized, see (15.9)-(15.10), such that c/6 is a 
parameter, fixed by the boundary conditions, especially A. The fixed points are Gaussian G, Random Periodic RP (the generalization of the RP 
fixed point for A = 0), Self- Avoiding Polymers SAP, and Unphysical U. 



The first line is indeed equivalent to (6.1) using A(n) = —R"{u). The second line is new and contains the terms 
induced by the KPZ term, i.e. the term proportional to A in (15.2). 

Equation (15.6) possesses the following remarkable property: A three parameter subspace of exponential func- 
tions forms an exactly invariant subspace. Even more strikingly, this is true to all orders in perturbation theory [5]! 
The subspace in question is (0 < ti < 1/ A) 



A(n) = — [a + he 



-\u^^^\u 



(15.7) 



The FRG-flow (15.6) closes in this subspace, leading to the simpler 3-dimensional flow: 



= a + 4a^ + 4ac + 46c 
dth = 6(1 + 6a + 6 + 5c) 
dic = c(l + 6a + 6 + 5c) 



(15.8) 
(15.9) 
(15.10) 



This flow has a couple of fixed points, given on figure 14. They describe different physical situations. The only 
globally attractive fixed point is SAP, describing self-avoiding polymers. This fixed point is not attainable from the 
physically relevant inital conditions, which lie (as fixed point RP) between the two separatrices given on figure 14. 
All other fixed points are cross-over fixed points. 

In the Cole-Hopf representation, it is easy to see why the exponential manifold is preserved to all orders. Let us 
insert (15.7) in (15.5). The complicated functional disorder takes a very simple polynomial form [5]. 



S=f Z{x,t){dt-V^) Z{x,t)- [ [ Z{x,t)Z{x,t'){aZ{x,t)Z{x,t') + bZ{x,tf + cZ{x,t'f) 

Jxt Jx Jt<t' 



(15.11) 



time 



Figure 15: The three vertices proportional to a, b and c in equation (15.1 1). 
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The vertices are plotted on figure 15. It is intriguing to interprete them as particle interaction (a) and as different 
branching processes (6 and c): Z destroys a particle and Z creates one. Vertex b can e.g. be interpreted as two particles 
coming together, annihilating one, except that the annihilated particle is created again in the future. However, if the 
annihilation process is strong enough, the reappearance of particles may not play a role, such that the interpretation 
as particle annihilation or equivalently directed percolation is indeed justified. 

One caveat is in order, namely that the fixed points described above, are all crossover fixed point, and nothing 
can a priori be said about the strong coupling regime. However this is the regime seen in numerical simulations, for 
which the conjecture about the equivalence to directed percolation has been proposed. Thus albeit intriguing, the 
above theory is only the starting point for a more complete understanding of anisotropic depinning. Probably, one 
needs another level of FRG, so as standard FRG is able to treat directed polymers, or equivalently the KPZ-equation 
in the absence of disorder. 

16 Perspectives 

More interesting problems have been treated by the above methods, and much more remains to be done. We have 
applied our techniques to the statics at 3-loop order [6] and to the random field problem [9]. An expansion in 1/N, 
(by now we have obtained the effective action [7]), should allow to finally describe such notorious problems as the 
strong-coupling phase of the Kardar-Parisi-Zhang equation. Many questions are still open. Some have already been 
raised in these notes, another is wether FRG can also be applied to spin-glasses. We have to leave this problem for 
future research and as a challenge for the reader to plunge deeper into the mysteries of functional renormalization. 
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